A unique non-Landau/Fermi liquid in 2d for high T c superconductivity 
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It is shown that the main features of the high T c phase diagram can be calculated as a function 
of doping in a simple, essentially unique non-Landau/Fcrmi liquid in 2d with quartic interactions. 
This depends on a single parameter < 7 < 1 which encodes the strength of the interaction at 
short distances. A new d-wave gap equation has solutions that fall under a superconducting dome, 
which terminates at the renormalization group fixed point. Optimal doping is estimated to occur 
just below 3/2-7T 2 . The scale for T c is set by the recently measured universal nodal Fermi velocity 
and lattice spacing, and is estimated to be 120A' < T c < 160A" for LaSrCuO. 



There remain many fundamental unanswered ques- 
tions in the theory of high T c superconductivity in the 
cuprates. How can one theory interpolate between anti- 
ferromagnetic (AF) and superconducting (SC) order, 
since AF requires repulsive interactions and SC attrac- 
tive? What is the precise mechanism that gives d-wave 
pairing? What is the nature of the pseudogap? The 
most important guide in tackling these problems is the 
non-Fermi liquid behavior, as emphasized early on by 
Anderson[l|. The reason is that if one insists on a local 
quantum field theory description in the rotationally in- 
variant long wavelength limit, relevant interactions which 
lead to a low energy fixed point are exceedingly con- 
strained and in fact rare. In this paper we describe an es- 
sentially unique non- Fermi liquid theory in 2d based on 4- 
fermion interactions which automatically has 50(5) sym- 
metry. It thus gives a microscopic model where the ideas 
of Zhang may be explored^. The requirements of a con- 
sistent local effective field theory for expansion around 
the Fermi surface and relevance of quartic interactions 
make this theory inevitable. This model was proposed 
earlier in some preliminary work by one of us0,3] which 
focused on AF order; at the time the d-wave SC proper- 
ties were not understood. The model is in line with ideas 
that emphasize the role of the Hubbard and Heisenberg 
models 5] . However our model cannot be simply derived 
by a naive scaling limit of the Hubbard model since the 
latter only has at most SO(A) symmetry^. Once the 
theory is formulated, the main properties of high T c fol- 
low naturally in an analysis involving only simple 1-loop 
calculations, which are accurate since at the fixed point 
the coupling is small sa |. Our model reveals that high 
T c SC may be a remarkably universal phenomenon that 
manages to realize some subtle theoretical loopholes, and 
its main properties follow from the existence of the low 
energy renormalization group (RG) fixed point. Such 
quantum critical points were emphasized by Vojta and 
SachdevQ. We begin by motivating our model with RG 
arguments and derivations of it in two limits. The rest of 
the paper explains how we calculated the phase diagram 
shown in Figure [TJ 

Consider free non-relativistic particles with energy 
e(k) that is rotationally invariant, e.g. e(k) = k 2 /2m*. 
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FIG. 1: Calculated phase diagram as a function of hole doping 
based on a single parameter < 7 < 1, set equal to 1 (in- 
finitely strong coupling at short distances). What is shown 
are solutions 5' s , 5' q of the AF and d-wave gap equations below, 
which are proportional to the critical temperature. T } 
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is sim- 
ply the RG scale. The AF transition point at Iiaf = 
first order. The SC transition at /i* = ^-r is second-order anc 
corresponds to the fixed point of the renormalization group 
hi = h(xi) S3 0.13 is not universal. 
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At finite density the Fermi surface is a circle in 2d as 
shown in Figure For a small band of energies near the 
Fermi surface, k = kj?(k) + p(k) as shown and the parti- 
cles have a linear dispersion relation e(k) ss ef ± Vj? |p|. 
Let a p correspond to particles and b p to holes. Then the 
effective hamiltonian is 

H = / (d 2 p) [(v F \p\ - £) a t flp + (v F \p\ + $)blb p ] 
«/|p|<Ac 

(1) 

where fx = ji — e F is zero at zero temperature. 

Our expansion around the Fermi surface is in the same 
spirit as in[j|, but we now depart from it by demanding 
a consistent local effective quantum field theory that re- 
produces the above H . Since the energy corresponds to 
massless particles with a linear dispersion relation, we 
identify an emergent Lorentz symmetry and describe H 
using a relativistic field theory. There are only 2 known 
candidates which differ in whether the lagrangian is first 
or second order in derivatives. The first order case re- 
quires a multi-component Dirac field, and since here a 4- 
fcrmion interaction is an irrelevant dimension 4 operator 
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FIG. 2: Expansion around a circular Fermi Surface. The 
diamond corresponds to a half-filled lattice. 

in 2d, it cannot lead to a non-Fermi liquid. Furthermore, 
Dirac fermions usually require additional properties of 
the Fermi surface (Dirac points) as in graphene, where 
the multi-components have to do with sub-lattices. 

The other possibility is second-order in both space and 
time derivatives with action 

S = Jdt d 2 x (d tX -d tX + - v 2 F V X ~ ■ V X + ) (2) 

The Fermi velocity plays the role of the speed of light 
which just serves to convert units of space and time so 
it can be set equal to 1. This form for a fermionic field 
is very unconventional, and to a particle physicist it ap- 
pears to violate the spin-statistics theorem. Since the 
above kinetic term is crucial to all that follows, we give 
the following compelling arguments in favor of it: (i) It 
correctly reproduces the desired effective hamiltonian ([I]) 
for particles and holes near the Fermi surface. (See be- 
low.) Thus the free theory is perfectly hermitian and 
unitary in momentum space, i.e. has no negative norm 
states, (ii) In the condensed matter context, spin is a fla- 
vor and for spin i particles we simply double the number 
of components xf ;Xf- Since there is a total of 4 fields, 
by Fermi statistics there is a unique 4-fermion interaction 
with hamiltonian density 

Hint = 8ir 2 g xf xfx+ X; ( 3 ) 

Repulsive interactions correspond to positive g. Since the 
field has classical scaling dimension ^ in 2d, the above 
operator has dimension 2 and is therefore relevant. At 
low energies the coupling flows to an interacting fixed 
point with non-Fermi liquid behavior, (iii) Although our 
model was originally motivated by expanding near a cir- 
cular Fermi surface, it can also be obtained from inter- 
acting itinerant lattice fermion models like the Hubbard 
model at half-filling, thus it can interpolate between the 
two Fermi surfaces shown in Figure[2] At strong coupling 
the latter is known to correspond to the Heisenberg anti- 
ferromagnet, with a low energy description in terms of 
an 0(3) non-linear sigma model for a field <fi constrained 



to have fixed length with lagrangian C = 8^$- 9^(?[9l[lf3|. 
In our model the </> order parameter is 4> = \~ ax + / 'V% 
(see below) and the constraint on cj) follows from the sim- 
ple constraint x~X + equal to a constant. Inserting this 
into the <f> action one obtains the second order action ([2]) 
up to some irrelevant operators. This was pointed out 
in 0] and explained in more detail below in connection 
with doping. There is no similar construction for Dirac 
fermions. 

There is no violation of the spin-statistics connection 
since spin is a flavor in this context and the Pauli- 
exclusion principle is built in from the fermionic nature 
of the x-fields. The issue rather has to do with unitarity. 
The mode expansion of the fields is 

X"(x,t) = J ^(ate-*- + 6pe*-*) 

x +(x,i) = y^(-^e-** + ape*«) (4) 

where uj p — \Zp^ and p ■ x = ui p t — p ■ x. The additional 
minus sign in the expansion of x + is chosen so that the 
canonical quantization relations of the fields leads to the 
usual canonical relations in momentum space for the a 
and &'s. The canonical hamiltonian of the theory is pre- 
cisely ((TJ. Introduce a unitary operator C that distin- 
guishes particles and holes: CaC — a, CbC = —b where 
C 2 = 1. Then x + = C(x~yC an d in terms of fields the 
hamiltonian is pseudo-hermitian: W = CHC\^. It was 
understood long ago by Pauli that a pseudo-hermitian 
hamiltonian gives a consistent quantum mechanics with 
a unitary time evolution and real eigenvalues. In the 
present context, pseudo-hermiticity has additional mean- 
ing with regard to the kinematics of the expansion around 
the Fermi surface since C distinguishes particles and 
holes. Conservation of the physical momentum k is only 
equivalent to conservation of p for processes where par- 
ticles are paired with particles and holes with holes. For 
the study of SC, these are of course the processes we are 
primarily concerned with. We are thus only interested in 
eigenstates which are also eigenstates of C , and for these 
H = fft. 

The SO(5) symmetry is easiest to see if one introduces 
an iV-component version with fields a = 1,..,N, 
which has Sp(2N) symmetry Q. For spin | particles, 
N = 2 and Sp(4) = 50(5). The 50(5) has an SU{2) 
spin subgroup and a U{\) charge that commutes with 
it. The ± indices on the fields x^ correspond to elec- 
tric charge. One can construct an SO(5) vector of order 
parameters $ = {<p x , <fi y , <fi z , (/)+ , <fi~) where <fi is an electri- 
cally neutral SU (2) vector and (f>f are Cooper pair fields 
of charge ±2 which are SU (2) spin singlets: 

<t> = x~<?x + /V2, 0^=x|x| 5 <K = Xj^Xf (5) 

The interaction can be written in a manifestly SO(5) 
invariant fashion: H m t = ® ■ 
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It is important to carry out the RG directly in 2c?. As 
usual, the RG prescription involves two energy scales, the 
cut-off A c and a lower running scale A. In many calcula- 
tions A is a lower cut-off. Since the coupling g has units 
of energy in 2d, we define <?(A) = A<?(A) where 3 is di- 
mensionless. The 1-loop beta function is — i^ta = 9~&9 2 
which has a low energy fixed point at g* — 1/8. To un- 
derstand the phase diagram as a function of doping, it is 
first convenient to introduce the variable x = X/g where 
the fixed point value is at 1, =8. We also introduce 
a variable xq that encodes the strength of the coupling 
at the cut-off: xq = l/<?o where g(A c ) = A c g . We as- 
sume that at short distances the coupling is strong, i.e. 
g > <?». It will also be useful to define 7 = (x* — xq)/x* 
which is a small parameter between and 1 The cou- 
pling at short distances can be arbitrarily strong, where 
infinite coupling corresponds to 7 = 1. Integrating the 
beta-function with this initial short-distance data gives 
a linear form that is specific to 2d 
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(6) 



and turns out to be important in connection with hole 
doping, which we now turn to. 

In the non-linear sigma-model description at half- 
filling, the order parameter cj) is constrained to have fixed 
length. As explained above this constraint follows from a 
constraint on the x fields: xf X] + Xjx'i = ihA c . Then 
one can show <p ■ <p = 3h 2 A 2 /2. Relaxing this constraint 
moves away from half-filling. Thus a measure of hole 
doping is the 1-point function h = —i(x~x + )/A c . The 
overall scale of ft as a doping variable can be justified 
by reintroducing the chemical potential and noting that 
near the Fermi surface it couples to the above operator. 
Including the 1-loop order g self-energy correction to the 
propagator and expressing everything in terms of the x 
variables one obtains 
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(7) 



In Figure Q] the straight line is a plot of A/A c in eq. © 
as a function of the above h for x — 0, i.e. 7=1. (For 
7 7^ 1 it is not exactly a straight line.) It crosses the h 
axis at /i» = 3/2tt 2 , which is the location of the RG fixed 
point. Below this line, the energy scale is such that elec- 
trons described by our x fields are strongly correlated and 
exhibit non-Fermi liquid behavior since the fixed point is 
not a free field theory. Only (x~X + ) 7^ and no symme- 
tries are broken at this line. The region below is what 
is normally called the pseudogap. The scale A can be 
associated with a pseudogap temperature T pg . 

The AF phase can be analyzed by a standard mean 
field analysis. One introduces an auxiliary field s for the 
order parameter ip and derives the effective potential for 
constant s by performing the functional integral over the 



X fields. Minimizing this effective potential with respect 
to s gives the gap equation 
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(8) 



For positive g, there are solutions due to the compen- 
sating minus signs. Since s has dimension 2, define 
s = 6 2 A 2 . Then <5 S is a solution to the equation 
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(9) 



When g is small enough, the solution flattens out with 
5 S w 1+. This behavior is unphysical since the gap should 
be zero when g is zero. The resolution of this puzzle in- 
volves regulating the infra-red divergence with the low 
energy cut-off A and interpreting the result with the RG. 
Setting s — 0, the gap equation can be approximately 
re-expressed as 1/ 17(A) = 8/A c . This shows that at 
g = A c /8, a consistent non-trivial solution is s = 0. We 
interpret this as a first-order transition where S s drops 
discontinuously to zero. In terms of x this occurs at 
x af — jq?- ■ Since S s is in units of the cut-off it is mean- 
ingful to rescale it and define S' s = -r-S s , where the scale 
factor is given in terms of x in ^ . In Figure [T] we show 
the solutions to the gap equation as a function of dop- 
ing h. The Neel temperature Tjv is proportional to the 
zero temperature gap as we will describe below. The AF 
transition occurs at Iiaf = h{xAF) = 3/4ir 2 when 7=1. 

Introducing constant auxiliary fields q for the SC or- 
der parameters <pf and repeating the above mean field 
analysis gives a gap equation for q 2 = q + q~ that has the 
same form as ([8]) with s 2 — * — q 2 . This leads to the ex- 
pected result: for repulsive interactions (positive g) there 
are no solutions, i.e. no s-wave SC. This shows that the 
AF and SC phases do not compete, and the SC phase is 
not simply related to the AF one by the SO(5) symme- 
try. When one goes beyond mean field and incorporates 
momentum dependent scattering of Cooper pairs near 
the Fermi surface, an attractive d-wave channel opens 
up. Introducing non-constant auxiliary pair fields, one 
can derive the momentum dependent gap equation 



duj d 2 k' 



(2tt) 3 



G(k,k') 



q(k') 



(lu 2 + k' 2 ) 2 + g(k') 2 



(10) 



where q + = q~ up to a phase. The kernel G is related 
to a particular 4-particle Green function specialized to 
Cooper pairs of opposite momenta ±k and ±k'. In a 
rotationally invariant theory one can expand in circular 
harmonics: 



G(k,k') = Y^G t {k,k') cast{6-G') 

e=o 

00 

q(k) = ^Tq e (k)cos£9 (11) 
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where k is the magnitude of k and 9 — 9' is the angle 
between k and k'. Performing a low energy momentum 
expansion at 1-loop one finds an attractive I = 2 channel: 
G 2 {k,k') = -8ir 2 g 2 k 2 k' 2 where g 2 = 4g 2 /25A 3 . The 
solution to the gap equation has the characteristic d-wave 
form 



g(k) = 5lk 2 cos 20 = S 2 Jk 2 x - kl) 



(12) 



where 8 q is a constant solution to the integral equation 
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(13) 



To plot the gap as a function of doping one needs to ex- 
press the constant g 2 in terms of x: g 2 = 25x -2(tl x / x ,y » 
and one sees that it changes sign at x — x*, which im- 
plies that the SC phase terminates at the RG fixed point. 
The other important feature of the above gap equation 
for S q is that SC turns off at a lower, non-universal value 
of Xi w (1 — 0.17)2:*. Thus the non-zero solutions of the 
d-wave gap equation fall under a dome x\ < x < x*. 
Comparison with experiments shows that this dome as 
calculated is too narrow; however there are many effects 
we have neglected that could broaden it, e.g. interplane 
coupling, disorder, or higher order corrections. Numeri- 
cal study of solutions of (|13[) indicates that optimal dop- 
ing occurs in the tight range between 0.13 and 0.15 where 
the upper bound is simply h*. Solutions to the d-wave 
gap equation 5' — A<5 9 /A c are shown in Figure [T] 

As in the BCS theory, one expects that the critical tem- 
perature T c is proportional to the zero temperature gap. 
Rather than developing a full-blown finite temperature 
formalism, one can study the effect of a small finite tem- 
perature by introducing a relativistic mass term in the la- 
grangian — m 2 x~X + where m = aT with a a dimension- 
less constant. The justification for this is that in the free 
theory this thermal perturbation reproduces the correct 
electronic specific heat Cy oc Tkp/vF- Requiring the 
overall coefficient to match the standard non-interacting 
result extended to 2d fixes a = it 5 / 4 /\/6 « 1.7. We 
have verified numerically that the solutions of the AF 
and d-wave gap equations with this additional mass term, 
lu 2 — > uj 2 + to 2 , vanish when to is too large, and the crit- 
ical value m c — aT c is indeed proportional to the zero 
temperature gap up to constants of order 1 that depend 
weakly on 7. Restoring fundamental constants and the 
Fermi velocity vf one then finds kgT c = ~ vphd' q A c . At 
optimal doping the constant c is in the range .5 < c < .7 
depending on 7 and 5' w 0.11 is largely insensitive to 
7. A similar formula applies to the Neel temperature 
Tjv with c ~ 1. The Fermi velocity vf in our rela- 
tivistic model is a fixed constant that plays the role of 
the speed of light which is universal in the sense that it 
is independent of the coupling and thus the doping h. 



Remarkably, at low energies a doping independent uni- 
versal nodal Fermi velocity has been observed in recent 
years (ll). This can actually be viewed as a prediction 
of our theory, where the nodal direction k from (0, 0) to 
(tt, 7r) shown in Figure [5] is what is important since this 
is the direction where the half-filling diamond is closest 
to the nearly circular Fermi surface just inside it. Using 
the results in[ll| we estimate vf ~ lAevA — 210km/s 
for LaSrCuO. This leads to an estimate of T c at opti- 
mal doping. The cut-off A c should be equated with the 
inverse lattice spacing a. A practical form is 



vf 

T c = c— ■ 650K 
a 



(14) 



where vf is in evA and the lattice spacing in A and we 
have set a to its estimated value of 1.7. With the above 
vp and a — 3.8A for the CuO square lattice, this gives 
T c in the range 120AT < T c < 160AT depending on 7. The 
maximum T c occurs around 7 = 1/4. The simple result 
(TTJJ provides some insights on how to possibly obtain 
a higher T c : shorten the lattice spacing or increase the 
Fermi velocity by possibly reducing the effective electron 
mass to*. 

Many more properties can be calculated, such as 
anomalous corrections to the specific heat, a resistivity 
linear in temperature, and electron correlation functions 
in the pseudgap region. These additional results and a 
more detailed exposition will be be presented in the com- 
prehensive article 12j . 
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